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A century ago, J. J. Thomson 1 showed that the 
scattering of low-intensity light by electrons was a 
linear process (i.e., the scattered light frequency 
was identical to that of the incident light) and 
that light's magnetic field played no role. To- 
day, with the recent invention of ultra-high-peak- 
power lasers 2 it is now possible to create a suffi- 
cient photon density to study Thomson scattering 
in the relativistic regime. With increasing light 
intensity, electrons quiver during the scattering 
process with increasing velocity, approaching the 
speed of light when the laser intensity approaches 
10 18 W/cm 2 . In this limit, the effect of light's 
magnetic field on electron motion should become 
comparable to that of its electric field, and the 
electron mass should increase because of the rel- 
ativistic correction. Consequently, electrons in 
such high fields are predicted to quiver nonlin- 
early, moving in figure-eight patterns, rather than 
in straight lines, and thus to radiate photons 
at harmonics of the frequency of the incident 
laser light , with each harmonic having its own 
unique angular distribution 5 7 . In this letter, we 
report the first ever direct experimental confir- 
mation of these predictions, a topic that has pre- 
viously been referred to as nonlinear Thomson 
scattering 7 . Extension of these results to coherent 
relativistic harmonic generation 1011 may eventu- 
ally lead to novel table-top x-ray sources. 

In this experiment, we used a laser system that pro- 
duces 400-fs-duration laser pulses at 1.053-/im wave- 
length with a maximum peak power of 4 TW. The 50- mm 
diameter laser beam was focused with an f/ 3.3 parabolic 
mirror onto the front edge of a supersonic helium gas jet. 
The focal spot is consisted of a 7- /mi FWHM Gaussian 
spot (containing 60 % of the total energy) and a large (> 
100 /im) dim spot. The helium gas was fully ionized by 
the foot of the laser pulse. A half-wave plate was used 
to rotate the axis of linear polarization of the laser beam 
in order to vary the azimuthal angle (</>) of observation. 
We define 9 = 0° as along the direction opposite to that 
of the laser propagation and 4> — 0° as along the axis 
of linear polarization. In a linearly polarized laser field, 
electrons move in a figure-eight trajectory lying in the 



plane defined by the axis of linear polarization and the 
direction of beam propagation. 

While the observation of harmonics in laser-plasma (or 
electron beam) interactions has been made by several 
groups 12-16 , that alone is insufficient to unambiguously 
identify nonlinear Thomson scattering and its underlying 
dynamics. Several other mechanisms might generate con- 
tinuum or harmonics under our experimental conditions, 
and, therefore, need to be isolated and discriminated 
from the signal generated by nonlinear Thomson scat- 
tering: (1) continuum generated from self-phase modula- 
tion of laser beam in gas, (2) harmonics generated from 
atomic nonlinear susceptibility of gas or, especially, from 
the ionization process 17 , (3) continuum generated from 
(a) (relativistic) self-phase modulation of laser pulse in 
the plasma, or from (b) electron-electron bremsstrahlung 
and electron-ion bremsstrahlung, and (4) harmonics gen- 
erated from the interaction of laser pulses with a trans- 
verse electron-density gradient 14 . 

The main focal spot of the laser pulse undergoes 
relativistic-ponderomotive self-channeling when high 
laser power and gas density are used 18 . Side imaging (6 
= 90°) of the 1st harmonic light (at the laser frequency) 
from nonlinear Thomson scattering shows that the laser 
channel has a diameter of <10 /jm FWHM. However, in- 
terferograms 18 show that the diameter of the plasma col- 
umn is about 100-200 /im, which is created by the wings 
with intensities > 10 15 W/cm 2 (the ionization threshold) . 
Therefore, the light generated from laser-gas interaction 
should be observed to originate from the entire region of 
plasma, rather than from the narrow laser channel. Re- 
sults of side imaging (6 = 90° and arbitrary <f) of the 
2nd and 3rd harmonics using a matching interference fil- 
ter (10 nm bandwidth) show that the signal is emitted 
only from the narrow laser channel. In addition, the im- 
ages of the harmonics have spatial distributions similar 
to the images of the 1st harmonic light, and their profiles 
vary in the same way as the laser power and gas den- 
sity are changed. This rules out the possibility that the 
harmonic signal observed in the side images is a result of 
laser-gas interaction ((1) and (2)). 

According to theory 5-7 , the harmonic signal generated 
from nonlinear Thomson scattering should have two im- 
portant features: (1) it is linearly proportional to the 
electron density because it is an incoherent single elec- 
tron process (the harmonics generated from a collec- 
tion of electrons interfere with each other destructively, 
leaving only an incoherent signal, which is equal to the 
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single-electron result multiplied by the total number of 
electrons which radiate), and (2) it increases roughly 
as 7™, where n is the harmonic number, and gradu- 
ally saturates when ao is on the order of unity 7 , where 
a = eE/m aj c = 8.5 x lO-^At/xnrj/^IW/cm 2 ] is the 
normalized vector potential, E is the amplitude of laser 
electric field, and I = cE 2 /8ir is the laser intensity. These 
are characteristically different from the behavior of any 
other mechanisms. For instance, bremsstrahlung radia- 
tion should be proportional to the square of gas density 
(N e ■ N e or N e ■ Ni). In this experiment, the intensity of 
the harmonic signal was determined from the peak inten- 
sity or the average intensity of the images of harmonics, 
when it was plotted as a function of the observing angle, 
gas density and laser power. Both showed the same vari- 
ations. Figure [l] shows the variation of the 2nd harmonic 
signal as a function of laser power and plasma (electron) 
density. The experimental results show a reasonable fit 
with the theoretical predictions. The 1st and 3rd har- 
monics show the same match with the theory. 

Although the above two observations are consistent 
with nonlinear Thomson scattering as the source of the 
harmonic signal, the observation of the unique angular 
patterns is necessary in order to prove that the detailed 
dynamics of nonlinear Thomson scattering are indeed the 
same as the theoretical prediction. Figure ||(a) shows the 
</> -dependence of the 2nd harmonic signal at 9 = 90°. 
The experimental results match qualitatively with the 
theoretical prediction, both having a quadrupole-type ra- 
diation pattern, which is characteristically different from 
the dipole pattern for other mechanisms (l)-(4), and lin- 
ear Thomson scattering. Other measurements such as 
the ^-dependence of the 2nd harmonic light at 8 = 51° 
(an "anti-dipole" pattern), shown in Fig.|^(b), and the 
^-dependence of the 3rd harmonic light at 9 = 90° (a 
"butterfly" pattern), shown in Fig.y, were also made, 
all showing reasonable matches between the experimen- 
tal data and the theoretical predictions. Such angular 
radiation patterns directly prove that electrons do in- 
deed oscillate with figure-eight trajectories in an intense 
(relativistic) laser field. The angular pattern of the 1st 
harmonic light (linear component )of nonlinear Thomson 
scattering is also included in Fig. 0(b) for comparison. 

Measurements of the spectra of the harmonics show 
that each of the spectra of 2nd and 3rd harmonics con- 
tains a peak at roughly the harmonic wavelength and a 
red-shifted broader peak, as shown in Fig.^J. The red- 
shifted broader peaks are believed to be part of the har- 
monic spectra generated by nonlinear Thomson scatter- 
ing, because they vary in amplitude proportionally with 
the corresponding unshifted harmonic signals when the 
gas density and the laser power are changed. It was ex- 
pected that the spectra of harmonics should be broad- 
ened tremendously for electrons in a high-fluid-velocity 
plasma wave 6 . A fast-phase- velocity electron plasma 
wave (with a maximum fluid velocity of as large as ~0.2 c, 



where c is the speed of light in vacuum) excited by stim- 
ulated Raman forward scattering 19 was observed in this 
experiment at high laser power and gas density. However, 
the fact that the spectral distribution of the harmonics 
was not observed to change significantly with variation 
of gas density and laser power, when the plasma wave 
amplitude was, indicates that such spectral structure has 
nothing to do with the collective drift motion of electrons 
in the plasma waves. Although the angular radiation pat- 
terns of the harmonics could also be affected by such a 0.2 
c fluid-velocity oscillation, the changes are not significant 
enough (compare the solid and dash lines in Fig.|2|(a)) 
to be identified from the experimental data 7 . In other 
words, all measurements done in this experiment match 
qualitatively with the prediction of incoherent nonlinear 
Thomson scattering of electrons without drift motion; 
the results appear not to be affected by the existence 
of plasma waves, probably due to destructive coherent 
interference. The absolute scattering efficiency is mea- 
sured to be 8 x 10~ 4 and 1 x 10~ 4 photons per electron 
per pulse for the 2nd and 3rd harmonics (including both 
the unshifted and red-shifted spectral components), re- 
spectively, at 8 = 90° , 4> — 50° , for an angle of collection 
of 7 x 10~ 3 steradians. These numbers match reason- 
ably well with the theoretical predictions for incoherent 
nonlinear Thomson scattering, which are 8 x 10 -4 and 
5 x 10~ 4 , respectively. 

In summary, the results reported here confirm for the 
first time several predictions of relativistic electrody- 
namic theory, which were formulated forty years ago, co- 
incident with the invention of the laser. As predicted 
7 , a century-old fundamental "constant," the Thomson 
cross-section, is now shown to depend on the strength of 
light. 
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FIG. 1. Intensity of the 2nd harmonic light at 6 — 90°, 
4> = 50° (a) as a function of laser pulse energy at 
6.2 x 10 19 -cm~ 3 electron density and (b) as a function of 
plasma electron density at 0.8-J laser pulse energy. (For 1-J 
laser pulse energy, the laser intensity is 4.4 x 10 18 W/cm and 
ao is 1.88.) Each data point represents the result of a single 
laser shot. The theoretical predictions for zero drift velocity 
are plotted in solid lines for comparison. The only fitting pa- 
rameter is just a constant for normalization in all figures. The 
laser pulse energy refers to the energy in the main focal spot. 
The inset shows the coordinate system used. 
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FIG. 2. Polar plots of the intensity of the 2nd harmonic 
light as a function of azimuthal angle (<j), in degrees) for 
0.8-J pulse energy and 6.2 x 10 19 -cm -3 electron density at 
(a) 9 = 90° and (b) = 51°. The intensity is in arbitrary 
units. The solid circles represent the experimental data. The 
solid and dash lines represent the theoretical results for zero 
and nonzero drift velocity, v = 0.2c, in the laser propagation 
direction, respectively. The open circles and triangles repre- 
sent the experimental data for the 1st harmonic signal taken 
at two different runs under the same conditions. The dotted 
line represents the theoretiocal result. Its dipole radiation 
pattern (peaked at (j> =90°) confirms that there is no depo- 
larization effect in the plasma and that the collective effect 
of plasmas on the angular pattern is not significant (at least 
outside of a narrow cone along the axis of laser propagation). 
Although the data for 4> — 180° ~ 360° are not plotted, it 
should be just a mirror image of the data for <f> = 0° ~ 180° 
because of the intrinsic symmetry of the laser field, as pre- 
dicted theoretically. Such symmetry has been verified in the 
experiment. 



FIG. 3. Polar plot of the intensity of the 3rd harmonic light 
as a function of azimuthal angle ((f), in degrees) at 6 — 90° 
for 0.8-J pulse energy and 6.2 x 10 19 -cm~ 3 electron density. 
The intensity is in arbitrary units. The solid circles represent 
the experimental data. The solid line represents the theo- 
retical result for zero drift velocity. The angular patterns of 
harmonics should not be sensitive to variation of laser inten- 
sity, as expected from theory and checked in the experiment. 
This is crucial to the success of our measurements because it 
alleviates the error caused by fluctuation of laser intensity. 
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FIG. 4. Spectra of the (a) 2nd and (b) 3rd harmonics at 
6 = 90°, 4> = 50° for 0.8-J pulse energy and 6.2 x 10 19 -c m - 3 
electron density. Vertical lines indicate the wavelengths of the 
unshifted 2nd and 3rd harmonics. The intensities are ploted 
in arbitrary units. The spectra do not change with variation 
of <f> at any specific 0, so the angular distributions measured 
are not affected by the choice of filter bandwidth. 
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